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Abstract. We have modelled the evolution of the number of Lya absorbers with redshift, resulting from the 
evolution of the ionising background and the Hubble expansion. The contribution of quasars (QSOs) and galaxies 
to the Hl-ionising UV background has been estimated. The QSOs emissivity is derived from recent fits of their 
luminosity function. The galaxy emissivity is computed using a stellar population synthesis model, with a star- 
formation history scaled on observations of faint galaxies at A > 1500A. We allow for three values of the fraction 
of ionising photons that can escape the interstellar medium, f^sc = 0.05, 0.1 and 0.4. The Intergalactic Medium is 
modelled as made of purely-absorbing clouds with the distribution in redshift and column density obtained from 
QSOs absorption lines. For the adopted values of /esc, the contribution of galaxies to the ionising UV background 
is comparable or greater than that of QSOs. Accounting for the contribution of clouds to the UV emission, all 
models with f^sc ^0.1 provide an ionising flux compatible with local and high-z determination, including those 
with a pure QSOs background. The observed z ~ 1 break in the evolution can be better explained by a dominant 
contribution from galaxies. We find that models in A-cosmology with nm=0.3, nA=0.7 describe the flat absorbers 
evolution for z ^ 1.0 better than models for f2m=1.0. 
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1. Introduction 

The evolution of the Lyman forest is governed by two main 
factors: the Hubble expansion and the the metagalactic 
UV background. At high redshift the expansion, which 
tends to increase the ionisation of the matter, and the UV 
background, increasing or non-decreasing with decreasing 
redshift, work in the same direction and cause a steep evo- 
lution of the number of lines with z. At low redshift, the 
UV background starts to decrease with decreasing red- 



dN/dz cx (1 + z)2.i9±o.27 1.5 < 2 < 4, while dN/dz oc 
(1 + z)" i6±o i6 at z < 1 (IWeymann et al. 1998| ). 



Recent numerical simulations have been remarkably 
successful in reproducing the observed evolution (see, for 
example Dave et al. 199£; Machacek et al. 2000| ), leav- 
ing little doubt about the general interpretation of the 
phenomenon. However the same simulations predict the 
break in the dN/dz power-law at a redshift z ~ 1.8 which 



shift, due to the reduced mimbei' of ionisiug sources, and 
this effect counteracts the Hubble expanaion. Ay a result 
the evolution of the number of lines slows down. 

In a recent study of the evolution of the Lyman forest 
Kim, Cristiani, fc D'Odorico (2001) have shown that the 



appears too high in the light of the new results of Kim 



et al. (2001). This suggests that the UV background im- 



number of Lya lines per unit redshift, dN/dz, is well de- 
scribed by a double power-law with a break at z 1. For 
column densities in the interval TVhi = 
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plemented in the simulations is not completely correct: 
QSOs have been considered as the main source of ionis- 
ing photons, and, since their space density drops below 
z 2, so does the UV background. However, galaxies can 
produce a non-negligible ionising flux too, perhaps more 
significant than p reviously assumed, as shown by recen t 
measurements by Bteidel, Pcttini, fc Adelberger (2001 ). 
The galaxy contribution can keep the UV background rel- 
atively high until at z ~ 1 the global star formation rate 
in the Universe quickly decreases, determining the change 
in the number density of lines. 
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In ihis paper we recompute the contribution of QSOs & Jakobsen 1990: Madau 1991, 1992) 



and galaxies to the UV background following the recipes of 



Madau (1991, 


1992); 


Haardt & Madau (1996 


) and 


Madau, 


Pozzet 


;i, & Dickinson (1998) (Sect. ^j). In Sect 


. 1^ the 



results are compared with constraints on the UV back- 
ground derived from the proximity effect and the Ha emis- 
sion of high galactic latitude clouds. In Sect. ^ the evolu- 
tion of the number of Lya lines per unit redshift is com- 
puted according to a simple analytical model for different 
relative contributions of galaxies and QSOs. In Sect. |^ we 
summarise the results and discuss some consequences of 
them in terms of present and future observations. 

All the computations are carried out for two flat cos- 
mologies: an Einstein-De Sitter cosmology {flm = 1, ^^A = 
0) and a A cosmology {ilm = 0.3, JIa = 0.7). We have 
adopted Hq^IQ km s^^ 



dFreedman et al. 200lD 



Mpc^i 

throughout and scaled to this value the data derived from 
the literature. 



2. The Ultraviolet Background 

The mean specific intensity of the Ultraviolet background 

J (^obsi ^obs 

), as seen at a frequency i^ohs by an observer at 
redshift ^obs, can be derived from 



J (^obs ) ^obs 



1 

47r 



(1 -I- Zobs) 

(1 + zf 



-e{v,z) 



dz 



(1) 



where v = i^obs(l + z)/(l + Zobs), e.{i',z) is the proper 
space-averaged volume emissivity, Tes (t'obs, ^obs, -2) is the 
effective optical depth at z^obs of the Intergalactic Medium 
(IGM) between redshifts Zobs and z, and dl/dz the prope r 
line element ( |Madau 199"l| , [l992| ; [Haardt fc Madau 1996| ). 

The emissivity e(i^, z) should include a contribution 
both from direct sources of UV radiation (e.g. QSOs and 
galaxies) and from the IGM clouds themselves, through 
continuum radiative recombination of the gas ( Haardt & 



Madau 1996). For the sake of simplicity, we consider here 
the catc ot a purely absorbmg IGM, thus omitting ra- 
diative recombmation. The ettect ot this omission will be 
discussed later. 

The line element can be written as 



Toff (A 

obs ; ^obs 1 ^ ) 



dz' 



dNu,fiNii„z'){l~e-^^^'y), (3) 



where f{N}ii,z') = d N/dNnidz' is the distribution of 
absorbers as a function of redshift and column density of 
the atomic hydrogen A'^hi, and r(A') is the optical depth of 
an individual cloud for ionising radiation at a wavelength 
A' = Aobs(l + 2obs)/(l + z). For 228A< A' <912A, the 
main contribution to absorption of UV photons is given 
by the ionisation of Hi, therefore 



t(A') = rHi(A') = A^HifTHi 



(4) 



V912A 

where cthi = 6.3 -10"^* cm~^ is the photo-ionisation cross- 



section at the Lyman limit for Hi (Osterbrock 1989). 

The ionisation of Hei at 504A is not considered: He I 
being almost completely ionised, its contribution to the 
total opacity is negligible (Haardt & Madau 1996). For 
completeness, we have included the contribution of He 11 
ionisation to the opacity, although it does not affect the re- 
sults presented in this paper. He 11 is ionised for A' <228A: 
thus, 



t(A') = Thi + A^HciiCTHo 



/ A' 



V228A 



(5) 



with cTHcii = 1-58 • 10"!^ cm-2 (Osterbrock 198£). The 



column density A^hc n can be derived from A^h i by solving 
the radiative transfer within a cloud. When clouds are 
optically thin at 228A, a simple solution can be found 
( [Haardt fc Madau 1996|) : 



AT, 



He II 



1.8 M 



J(912A) 
J(228A) ■ 



(6) 



While Eqn. (^) does not hold for optically thick clouds, its 
use in the optically thick case can still provide correct esti- 
mates for the cosmic opacity, when Aobs > 228A ( Haardt] 



dl 

dz 



H{z){i + zy 



(2) 



fc Madau 1996). Since the He 11 contribution to Toff de- 
pends on the UV background, we have solved Eqn. (|l]) 
iteratively for every value of Zobs- 

For the redshift and column density distribution of ab- 
sorption lines we have adopted the usual form 



where c is the velocity of light and H{z) — iJo[^^m(l + 
z)^ + JIa]^/^ is the Hubble parameter for a flat universe 

{n^nrn + nA - 1). 



f{NK„z) 



A^Hi 



A 



1017 / V 1017 cm- 



a + zy 



(7) 



Fits of the absorption line distribution show that the in- 



2.1. Opacity 

The effective optical depth Toff through the IGM is defined 
as e"'^"" = (e^^), where the mean is taken over all the lines 
of sight from the redshift of i nterest. For a Poisson distri 
bution of discrete absorbers ( Paresce et al. 1980 ; MoUer 



dex (3 varies for different ranges of A^hi (Fardal et al 
1998; Kim et al. 2001). However, it is possible to de- 
scribe the cloud distribution with /? — 1.46 over several 



decades in A^hi ( Petitjean et al. 1993 ). For simplicity, 
we have adopted this single val ue over the whole c olumn 
density range considered here. Kim et al. (2001 ) com- 



bined high resolution VLT/UVES observations of 3 QSOs 
with literature data and derived a line number density 
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A 


/3 


7 


iVHi/cm-2 






0.50 


1.46 


0.16 


10" - 1.58 


10^' 


z < 1 


0.13 


1.46 


2.19 


10" - 1.58 


10" 


z > 1 


0.17 


1.46 


1.55 


1.58 • 10^^ - 


10^0 





Table 1. Parameters adopted for the distribution of ab- 
sorbing clouds (Eqn. |^). 



per unit redshift dN/dz = 9.06(1 + z)^ !^, for the col- 
umn density range iVni = lO^'^'^'^^^^cm^^ and z > 1.5. At 
lower redshifts. HST observations from the QSO absorp- 



tion li me key proiect show a slower evolution with red- 



effective optical depth Toff (Aobs, -^obs, = 1, as a function 
of Zobs- For radiation at Aobs = 912 A Az decreases from 
1.8 at Zobs=0 to 0.08 at Zobs=5. The same trend can be 
seen for Aobs = 600A but with larger values of Az (less 
absorption), because of the dependence of the Hi ionisa- 
tion cross-section on the wavelength. Since only radiation 
from local sources can easily reach Zobs, it is not necessary 
to compute the integral in Eqn. (^) up to z = cjo (or to the 
maximum z for which UV emitting sources are available). 
We have used in our calculation z,„ax=5- 

The absorber distribution we have adopted produces 
opacit ie s that are very s imila r to those of Haardt fc Madai; 
^19961) 



Fardal et al. 



(1998 



also 



ShuU et al. 19991) 



shift, with dN/dz = 34.7 (1 + z)"-^^ dWeymann et al 
19 98|)~|The chan ge in evolution, according to the results 



of Weymann et al. (1998 ) and Kim et al. (2001 ), occurs 
at z ~ 1. In this paper, we use 7=2.19 to describe the 
redshift evolution of the Lya forest (lO^'^ < A^Hi/cni^^ < 
1.58 • 10") at z > 1, and 7=0.16 at z < 1. By choos- 
ing A = 0.13 and A = 0.50, the integral of Eqn. (^ 
over A^Hi = 10^^'^''~^^cm^^ reproduces the results of Kim 



et al. '2001) and Weymann et al. (1998), respectively. 
The distribution of Lyman Limit systems is derived in an 



analogous way from Storrie-Lombardi et al. (1994). For 
^7 



1.58 • 10^^ < A^Hi/cm-^ < 10^", we use A = 0.17 and 
7=1.55. The parameters adopted for the distribution of 
absorbers are summarized in Tab. |^ 



As already pointed out by many authors (Madau 1991 



have derived values for A, (3 and 7 by fitting the dis- 
tribution of absorption lines in several ranges of column 
density. The opacity provided by their model is smaller 
than the one presented here. For example, at Zobs=3 we 
reach Tcff (912A) = 1 for Az = 0.18 while it is Az = 0.25 
for model A2 in |Fardal et al. (1998 ) . For the same emis- 
sivities, the opacity of Fardal et al. (1998) will result in a 
UV background higher than ours. The difference increases 
with Zobs and reaches 0.1 dex at Zobs=3. 

2.2. QSO emissivity 

The QSO contribution to the UV emissivity has been de- 
rived from the QSO luminosity function, for which we have 
adopted the double power-law Pure Luminosity Evolution 



1992; Haardt fc Madau 1996), the UV background be- model (Boyle et al. 1988) 



comes more dominated by local sources as the redshift 
increases. This can be seen in Fig. |l|, where we show the 
distance in redshift Az — z — Zobs corresponding to an 



0(L,z) 



L*(z) 



LHz) 



(8) 




where /3i and are the faint- and bright-end of the lumi- 
nosity function. A few functional forms have been adopted 
for the redshift evolution of the break luminosity L*{z). 
Using a sample of over 6000 QSOs with 0.35 < z < 2.3, 
Boyle et al. (2000) find that the evolution is well fitted by 
a second-order polynomial of the form 



L*(z) = L*(0) 10 



k\z-\-k2Z 



(9) 



0.1 - 



Fig. 1. Distance in redshift Az between an observer 
at Zobs and a point at z = Zobs + Az for which 
Teff (Aobs, Zobs, z) = 1, as a function of Zobs- Plots are 
shown for radiation observed at the Lyman edge and for 
Aobs = 6OOA. The plot does not include the contribution 
to the opacity of the He 11 ionisation, that will affect the re- 
sults at 6O0A. Assuming J(912A)/ J(228A)=50 (the mini- 
mum value in our simulations, when only QSOs contribute 
to the background), Az at 6OOA will decrease to « 2 in 
the range < z < 0.5. 



At redshifts z ~ 2, the luminosity evolution stops and the 
comoving number density of QSOs remains constant up 
to z ~ 3 ( iBoyle et al. 1988| ). At z > 3 the QSO num- 
ber density declines dramatically. A recent study by Fan 
|et al. (20"01 ) on a sample of 39 high-redshift QSOs from the 
Sloan Digital Sky Survey, suggests that the number den- 
sity declines as e~^''^^^ in the redshift range 3.6 < z < 5. 

Boyle et al. (2000) provide the best-fitting parameters 
of the B-band luminosity function for the two cosmologies 
adopted in this paper. The B-band proper emissivity can 
then be derived through the integral 



;(z.B,z) = {l + zf j^^^,^L(j){L,z)dL 
= {\ + zf 10'=i^+'=22=' e{vB,Q): 



(10) 
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where the factor (1 + z)^ is used to transform comov- 
ing into proper number densities. For the flat Einstcin- 
De Sitter universe and the A-cosmology, we have derived 
e{vB,0) = 1.2-1024 and 7.1-1023 /lyg erg s'^ Hz"! Mpc'^, 
respectively. Lmin scales with z as L* and is chosen to cor- 
respond to M^'^^ = —22 at z = 3, for both cosmologies 
(M™'''^ « -18 at z = 0; |Madau 1992] ). Since Eqn. (|) 
well describes a flat luminosity evolution for z ~ 2, we 
have adopted the emissivity of Eqn. (10) up to z = 3. 
For z > 3 we have adopted the exponential decline of 



Conversion factors for other choices of the IMF and stellar 



libraries among those provided by the Bruzual & Chariot 
model typically differ over a range of w 0.3 dex. The same 
span is also observed when comparing conversion factors 



obtained with different models (Kennicutt 199S). 

The star-formation history at z < 2 has been cal- 
ibrated on the UV comoving emissivities tabulated by 
Madau et al. (1998) for a flat rim=l universe. Emissivities 
at A = 2800A have been derived from the Canada-France 



et al. ( pOOlD . 
Boyle et al. 



Fan 



(2000) derived the R-ba,nd luminosity func- 



Rcdshift Survey (Lilly et al. 1996) in the range 0.2 < z < 1 
and from a HDF-north sample of obje cts with optical 
photometric redshifts for 1 < z < 2 ( Connolly et al 



tion hiim observations in the QSOs UV restframe. apply- 
ing the K-correction f or the composite QSO spectrum of 



Cristiani fc Vio (199C ) . For consistency, we have used the 



same spectrum to derive the UV emissivity for A > 1050A. 
For A < 1050A we have used a power-law, e{v) oc v"^-^, as 
meas ured on a sample o f radio-quiet QSOs observed with 
HST dZheng et al. 1997|). 



Wei ha^re also derived the QSO cm issivitjr from the 
work of La Franca fc Cristiani (1997 ). They fitted the 
luminosity function on a sample of 326 objects (the 
Homogeneous Bright QSO Survey) by using a different 
luminosity evolution. Results obtained with this emissiv- 
ity are very similar to those for the emissivity discussed 
above and are not presented in this paper. 



2.3. Galaxy emissivity 

The galactic emissivity in the ionising UV has been de- 
rived following the method outlined by Madau et al 
|(1998| ). The comoving UV emissivity at A > 1500 A (rest 



1997| ). [Madau et al. (1996| ) have derived the emissivity at 
A = 1500A from objects selected in two redshift ranges at 
z '--^ 3 and z ~ 4, from object selected on the HDF-north 
with the UV dropout technique. The lower emissivities at 
z > 2 suggested that the star formation rate has reached 
a maximum at z ^ 1 — 2. However, a different picture 
emerged from the ground-based survey of ^teidel et al. 
(1999|) . They selected Lyman break galaxies on a larger 



frame) |can be derived trom galaxy surveys as a tunction 
of the rcdshitt. iJecausc U V light is mamly produced by 



area than the HDF, refining the colour selection criteria 
with spectroscopic observations of a few object in the sam- 
ple. The derived emissivity (at A — 1700A) for z ~ 3 is 
still consistent with the Madau et al. (1996| ) value, while 
the emissivity at z ~ 4 does not show a steep decline. A 
value for the emissivity at A = 1700 A in the redshift bin 
2.5 < 2 < 3.5 can be derived from the luminosity function 
fitted by |Poh et al. (200l| ) on a combined ground-based 
and HST database. The datapoints are shown in Fig. |^ as 
a function of redshift. We have corrected the emissivities 
of [Steidel et al. (1999D and |Poli et al. (200 1| ) to include all 
objects with luminosit ies from to oo. Th e larg e errors in 
the da ta derived from Steidel et al. (199S| ) and Poll et al 



(2001 ) reflect the unc ertainties in the faint -end slope a 
of the Schechter (1976 ) luminosity fun ction, ^tcidcl et al 



sivity into a star-tormation history ot the universe. It we 
assume that the mean luminosity evolution of the galax- 
ies in the universe can be described with a single galactic 
spectrum compatible with the derived star-formation his- 
tory, a stellar population synthesis model can be used to 
derive the emissivity at any wavelength. 



short-l |ved UU stars, it is possible to convert the U V cmis- (199^ derive a = -1.6 ± 0.13 and |Poh et al. (200l|) 



We have used the latest version of the Bruzual & 



Chariot (1993 ) stellar population synthesis models, up- 



dated with a new set of stellar evolutionary tracks and 



a = -1.37 ± 0.20. In the [Madau et al. (1998]) tabulation, 
a — —1.3 was used. 

A smooth star-formation history has been derived from 
the observed UV emissivities in Fig. |[ by using Eqn. (11) 
and correcting for dust internal extinction according to 
the Calzetti 's ( 1997[ ) attenuation law. At high redshift, we 
have adopte d the flat star-forma tion history suggested by 
the work of Steidel et al. (1999). Synthetic galactic spec- 
tra have then been produced with the Bruzual & Chariot 



spectral libraries (Bruzual fc Chariot 2001 ; see also Liu (200T| ) code. In Fig. || we show the evolution of the mod 



et al. 2000[) . A [Salpeter (1955[ ) IMF with 0.1 < M^/Mq < 
100 has been adopted. We have chosen the quasi-empirical 
library of stellar spectra for solar metallicity, derived from 
observations for A > 1150A and from model of stellar 
atmospheres for A < 1150A. The conversion factor be- 
tween UV luminosity and SFR has been derived from the 
adopted stellar model, in the limit of continuous star- 
formation (Madau et al. 1998; Kennicutt 1998). For ra- 
diation at A — 1500 A we obtain 



SFR 



o7 erg s 
= 7.9 - 10^^ -4— 



Hz- 



M 



yr" 



(11) 



elled UV emissivity at 1500 A (solid fine) and at 2800A 
(dotted fine), for a colour excess E(B-V)=0.1. The model 
is compatible with both data at 2800A and z < 2 and 
data at 1500A (and 1700A) and z > 2. For the chosen 
E(B-V), the model is also consistent with the evolution 
of the emissivity in the optical-NIR regime for z < 2, as 
tabulated by Madau et al. (1995, not shown in Fig. |^). 

Unfortunately, it has not been possible to repeat 
the same procedure to model the emis sivity (and star- 
formation history) in A-cosmology. Only Poll et al. (2001 ) 
present a lu minosity function de rived assuming Qm = 
0.3, JIa — 0.7. Steidel et al. (1999 ) give the emissivity for 
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Fig. 2.Comoving UV emissivity in a 
flat Qm = 1 universe, from literature 
(dots) and from the spectral synthesis 
model described in Sect. 2.3 (lines). The 
model emissivity has been corrected 
for inter na l dus t extinction, using the 
Calzetti's (1997) attenuation law with 
E{B -V) = 0.1. Because of Eqn. (0) 
and of the redshift-independent dust 
correction, emissivity at A = 1500A 
is proportional to the SFR history 
adopted in the model (scale on right or- 
dinate). 



the objects visible in their survey, but they do not provide 
a luminosity function for an extrapolation to fainter lumi- 
nosities. Therefore, we have used the model of the emis- 
sivity for r2m=l and scaled it with a redshift-dependent 
correction: for sufficiently small redshift bin, it can be 
shown that the ratio of the emissivities in the flat A- and 
Einstein-De Sitter cosmologies is ^^0. 7 -f 0.3(1 -I- z)^/{l + 
z)^'^. The e missivity derived in this way is consistent with 
the data of |Poh et al. (200l| ). Because of Eqn. (|ll|), the 



Sest /esc « 5% ( 


Leitherer et al. 1995; 


Hurwitz et al. 1997 


Beckman et al. 2001 


). Steidel et al. (2001 


) analised a com- 



same ratios applies to the star-formation histories. 

The synthetic spectrum has then been used to derive 
the emissivity for the ionising UV. Due to the absence of 
observations, synthetic s pectra rely on models of stel lar at- 
mospheres for A < 9 12 A. Chariot fc Longhetti (2001 ) com- 
pared stellar spectra from different models and concluded 
that discrepancies on the ionising flux are not higher than 
0.1 dex. The uncertainty on our modelled emissivity also 
depends on the uncertainties in the determination of the 
SFR. To quantify the uncertainties in the adopted model, 
we have computed the effect on the emissivity of the vari- 



ation of the basic ingredients of the Bruzual & Chariot 
(200l|)]model (IMF, metaUicity, stellar libraries). We have 
used the same description for the UV emissivity at 1500A 
as in the main model (solid line in Fig. |^) and converted it 
into a star-formation history by using a conversion factor 
appropriate for the selected IMF and stellar spectra. The 
synthetic spectra obtained in this way typically differ by 
less than 0.2 dex at the ionisation limit and 0.3 dex at 
600A. 

Spectra at A < 912A also need to be corrected for 
the internal absorption by the galaxy interstellar medium. 
We describe this correction with the parameter /esc, i-e. 
the fraction of Lyman-continuum photons that can escape 
into the IGM without being absorbed by the interstellar 
medium, either gas or dust. A wide range of values can be 
found in the literature for /esc, derived both from models 
of radiative transfer and observations of Hi recombina- 
tion fi nes (5% < /esc < 60%; for a review, see Barkana 
& Loe 3 2001 ). UV observations of local starbursts sug- 



posite spectrum of 29 Lyman-break galaxies at z ~ 3.4. 
They derived a ratio between the flux densities at 1500 A 
and 900A /1500//900 = 4.6 ± 1, after correcting for the 
differential absorption due to the intervening IGM. The 
/1500//900 ratio for the unattenuated synthetic spectrum 
that we have used is very similar, /1500//900 ~ 5.3. If we 
assume that 40% of the radiation at 1500A is absorbed 



by dust (as obtained from the Calzetti's attenuation law 
with E(B-V)=0.1), the observed /1500//900 ratio is equiva- 
lent to /esc ~ 0.4 (if the internal absorption in the Lyman 
continuum do not change significantly with A). Because 
of the increase of the disk density with the redshift, /esc 
is expected to decrease with z; it is also found to depend 
heavily on the details of the distribution of the sources 
and the gas, i.e. whether the stars and/or gas are clumped 
or not (Wood & Loeb 200C). In this work, we will use a 
wavelength and redshift independent /esc, by which we 
multiply the synthetic spectrum at A < 912A. We will 
show results for /esc — 0.05 and 0.40, to cover the range 
of values suggested by local and z ^ 3.4 observations, and 
for an intermediate value, /esc = 0.1. 

Finally, the galactic emissivity has been converted 
from comoving to proper, multiplying by (1 -f zj^ . The 
total emissivity in Eqn. (^ is the sum of the QSOs and 
the galaxy contribution. 



3. The ionising background at 912A 

In Fig. H we show the modelled UV background, J{i/,z), 
at the Lyman limit as a function of redshift (solid lines) 
for the flat universe with r2m=l. The total background 
is shown as the sum of the QSO contribution (the same 
in each model; dotted fine) and the galaxy contribution 
(scaled with /esc; dashed fines). As predicted by other 
authors, for large values of /esc the ionising background 
produced by galaxies dominates over the flux from QSOs 
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Fig. 3.UV background at A = 912A for 
the models with /esc =0.05, 0.1 and 0.4 
(soUd lines), in a flat fim = 1 universe. 
Also shown are the separate contribu- 
tion of the QSOs (dotted line) and of 
the galaxies (dashed lines, each corre- 
sponding to a value of /esc)- The shaded 
area refer to the Lyman limit UV back- 
ground estimated from the pr o ximity 



effect ( Giallongo et al. 1996 ; []ook6 



et al. 1997; Scott et al. 2000). The ar- 



row shows an upper limit for the loca l 
ionising background (Vogel et al. 1995). 
The datapoint at z = 3 is derived from 
a compo site spectrum of L yman-break 
galaxies (Steidel et al. 2001 ). The mod- 



els and the datapoint of Bteidel et al 



(2001) have been multiplied by a z- 
dependent factor, to take into account 
the cloud emission (see text). 



( piallongo et al. 1997t pevriendt et al. 199S| ; |Shull et al 
199g|)7] 

At high redshift, the value of the UV background 
is constrained by the analysis of the proximity effect, 
i.e. the decrease in the number of intervening absorp- 
tion lines that is observed in a QSO spectrum when ap- 
proaching the QSO's redshift (Bajtlik et al. 1988). Using 
high resolution spectra, Giallongo et al. (1996| ) derived 
J(912A) = 5.0li '^ ■ IQ- ^^ erg cni'^ s'^ Hz'^ sr-^ for 
1.7 < z < 4.1 (see also Giallongo et al. 19991). Large r 
values are obtained by Cooke, Espey, & Carswell (1997), 



J(912A) = l.Otoj • 10"^^ erg cm"^ s^^ Hz^i sr^^ for 
2.0 < z < 4.5. A recent re-analysis of moderate resolu- 
tion spectra by ^cott et al. (200C| ) has lead to J(912A) 



7.0 



-1-3.4 



10 



-22 ^1 — 

erg cm s 



Hz sr for the same red- 



shift range. We show these measurements in Fig. y with 
a shaded area: the spread of the measurements obtained 
with different methods and data gives an idea of the uncer- 
tainties associated with the study of the proximity effect. 
Measurements of the ionising background at low red- 



shift are not less uncertain. Kulkarni & Fall (1993) re- 
ported a first tentative detection of the proximity effect in 
a sample of 13 QSOs at z < 1 observed with HST. They 



obtained J(912A) = 6^4' • lO"^"* erg cm"^ s"^ Hz" 



have shown that radiative recombination provides an im- 
portant contribution to the ionising background. Using 
only QSOs as source of ionising radiation, they obtain 
J(912A) = 5 • 10^22 ergcm-2 s'^ Hz"! sr"i at z=2.5, 
well within the shaded area in Fig. ^. Similar results are 
obtained by Fardal et al. (199"^ ). Although the cloud con- 
tribution to the background depends on the adopted emis- 
sivity and intergalactic absorption, we have obtained a 
rough estimate of its importance by using UV background 
sp ectra kindly provided by F. Haardt. At z=3, the models 
of iHaardt fc Madau (199"^ ) are a factor 1.7 higher than for 
the case of a purely absorbing medium; at z=0, the fac- 
tor reduces to 1.3. A simple linear interpolation between 
these two points reproduces the actual data for < 2; < 5 
within 5%. The models shown in Fig. ^jare multiplied for 
this z-dependent factor. 

Several models can produce a value of J(912A) com- 
patible with one of the measurements from the proxim- 
ity effect at z ~ 3 shown in Fig. |3[ from a a simple 
QSO-dominated background to models with f^sc ~ 0.2. 
However, a more stringent condition, fcsc ^ 0.10, is re- 



quired not to exceed the local upper limit of Vogel et al 



(1995). Bteidel et al. (2001) used their composite spectrum 



the large uncertainties due to the small number of avail- of Lyman-break galaxies and the UV emissivities of Steidel 



able alsorbers. Vogel et al. (1995|) derived a 2a upper limit 
J(912A)< 8.0 -10-^^ ergcm"^ s"^ Hz"^ sr~\ by study- 
ing Ha emission in a high latitude Galactic cloud. This 
upper limit is shown in Fig. |[ 

We must remember here that our model does not 
take into account the Lyman-continuum emission from 
recombination in Lya clouds. Haardt fc Madau (1996 ) 



et al. (1999[ ) to derive the ionising flux at z = 3. They ob- 
tained J(912A) = 1.2±0.3 10^21 erg cm'^ s'^ Hz^^ sr^^, 
a value consistent with our model for f^sc — 0.4. If Lyman- 
break galaxies with a spectrum similar to that observed by 
^teidel et al. (20"01 ) dominate the UV background, they 
will produce an ionising flux higher than the local and 
high-redshift estimates. 
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10 



13.64-16 



cm ^, for the two cosmologies adopted in 



this paper. Dotted lines refer to the evolution compatible with an ionising UV background due only to QSOs. Solid 
lines show the evolution when both QSOs and galaxies contribute to the background, for the models with /esc=0.05 
(upper line), 0.1 and 0.4 (lower line). Data points come from several observations in the literature for the column 
density range A'hi = lO^^-^"*"^^ cm~^, as given by Kim ct al. (2001 ). The modelled evolution has been normalized to 



the observed evolution in the redshift range 2 < z < 3. 



Similar results are obtained in the rim=0.3, f2A = 0.7 
universe (not shown). Because of Eqn. (||) and of the factor 
we have used in Sect. 2.3 to derive the emissivity for the 



A-cosmology, the contribution of galaxies to the emissivity 
is exactly the same in both the universe models adopted 
here. The QSOs contribution, instead, depends on an in- 
dependent fit of the luminosity function (Sect. ^.2[ ). In the 
new cosmology, the UV background produced by QSOs 
is slightly larger than the value for the Einstein-De Sitter 
model, by about 30% at the peak of the QSOs contribu- 
tion. [] 

4. The evolution of the Lyman Forest 

Further constraints can be obtained by studying the ef- 
fect of the evolution of the ionising back ground on the 
evolution of absorbers. Dave ct al. (1999 ) have studied 
the evolution of the low-redshift Lya forest in a hydrody- 
namic cosmological simulation, adopting an UV ionising 
background with the same redshift evolution as that of 



an analytical approximation for the evolution of dN / dz 
with Jy and the Hubble expansion. For clouds in photo- 
ionisation equilibrium with the background, it is easy to 
show that the evolution of lines above a given threshold 
in column density can be written as ( Dave et al. 1999| ) 



dN 
dz 



JVh,>JVS 



1 

Wz) 



rHi(2) 



(12) 



where H{z) is the Hubble parameter and (3 the coefficient 
of the power-law distribution of clouds with column den- 
sity (Eqn. 1^). T-i{i[z) is the photo-ionisation rate 

47r J(zy, z) 



hv 



aui{v)dv, 



(13) 



with i/Q the frequency of the Lyman limit and a]ii{v) — 
'^'Riivol vj^ the Hi photo-ionisation cross-section. 

In Fig. ^ we show dN/dz for Lyman forest clouds in the 



column density range iVui = 10 



13.64-16 



cm . Data points 



come from several sources in the literature and from new 



high resolution VLT/UVES spectra of three QSO (Kim 



at z = 1.7 for the number density evolution, dN/d'z. The 
change in evolution is primarily due to the drop in the 
UV ionising background, resulting from the decline in the 
QSO population. The formation of structure by gravita- 
tional growth plays only a minor role in the evolution. In 
the absence of structure evolution, is it possible to derive 



Haardi| Madau (1^^6r^ They found a sharp UansiLion 200l|). For each of our models, we have computed 



^ As pointed out by the referee, the background is, in prin- 
ciple, independent on the adopted cosmology. The slight dif- 
ference between the two QSO backgrounds arises because of 
the cosmology-dependent correction for incompleteness in the 
luminosity function. 



rHi(2;) and we have derived the evolution of dN/dz accord- 
ing to Eqn. (|l2|), for the two cosmologies adopted in this 
paper. The evolution has been normalized to the observed 
values for 2 < z < 3. In this redshift range, the UV back- 
ground is nearly flat for any of the models and it is easy to 
show, from Eqn. that dN/dz (x (1 + zf^'^-^, inde- 

pendently of the cosmology. By fitting the observed data, 
|Kim et al. (200l| ) have derived dN/dz oc (1 -t zf-^'^ for 
z > 1.5. A value /3 w 1.7 can well reproduce the evolution 
derived from the observations. This is consistent with fits 
of the density distribution, that give fi — 1.68 ± 0.15 over 
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Hi 



at z - 1 (Kim ct al. 2001). On the 



other hand, weaker Hnes (iVni ^ lO^cni are known to 
have a flatter distribution in column density (/3 1.4-1.5; 



not exceeding the upper limit for the local background at 
z = 0. 



Giallongo et al. 1996; Kim et al. 2001). This will produce a 



slower rcdshift evolution, as observed in this column den- 
sity range for z > 1.5 (7 
here that in Sect 



1; 



2.1 



Kim et al. 2001). We remind 



we have used /3 — 1.46, that provides 
a good description of the density distribution over a much 
larger column density range. 



The analysis of Kim et al. (2001 ) shows that the change 
in evolution occurs a t z sa 1, rather than at z ~ 1.7, as 
previously suggested ( Weymann et al. 1998 ). In Fig. ^the 
break at z « 1 can be reproduced if the contribution of 
galaxies to the background is dominant. This is because 
of the rapid decrease of the star-formation rate (and of 
rHi(z)) below this redshift (Fig. The photo-ionisation 
rate of a pure QSOs background, instead, peaks at z 2.5 
and has a slower evolution with z. It is interesting to note 
that the modelled evolution is closer to the observed data 
for the A-cosmology. For the Einstein-De Sitter universe, 
dN/dz grows for z < 1, which is not observed. However, 
the discrepancy may be mitigated when the effect of the 
formation of structures on dN/dz is taken into account 
( [Dave et al. 19991 ). 

The modelled dN/ dz do not depend on our approxima- 
tion of a purely absorbing intergalactic medium, since the 
contribution of cloud emission to Fh/z) is nearly consta,nt 



with rfldshift (see Fig, fi in Haardt Madau 199f ^ 



The analysis of Kim et al. (2001 ) shows that the break 
in the evolution of the Lya forest, dN/dz, occurs at z ~ 1. 
We used the evolution of the ionising background in our 
model to derive dN/dz, assuming that the formation of 
structure plays only a minor role ( Dave et al. 1999| ). The 
rapid decrease of star-formation for z ^ 1 can easily ex- 
plain the observed break, while a QSO-only background 
would produce a break sA z ^ 2 (Dave et al. 199!; ). For the 
ionising background to have an evolution similar to that 
of a galaxy-only model, high values of /esc are needed. 
This result apparently pushes the galaxy contribution in 
an opposite direction with respect to the estimates of the 
absolute value of the ionising background. The flat evolu- 
tion of dN/dz at z < 1 is much better reproduced adopt- 
ing a A-cosmology (r2m=0.3,r2A=0.7) rather than a rim=l 
Einstein-De Sitter universe. However, such a result needs 
to be confirmed after properly taking the formation of 
structure into account. 

A significant contribution from galaxies to the ion- 
ising metagalactic flux would correspond to a softening 
of its spectrum with respect to a purely QSO-dominated 
background. This effect would be particularly important 
at very high and very low redshift. Observations of the 
evolution of the Si iv to C iv ratio with redshift ( 3avaglic 
et al. 1997; Songaila 199^ ) seem to confirm the progressive 



5. Conclusions 

In this work we have derived the H i-ionising background, 
resulting from the integrated contribution of QSOs and 
galaxies, taking into account the opacity of the interven- 
ing IGM. We have modelled the IGM with pure-absorbing 
clouds, with a distribution in column density of neutral 
hydrogen, A'hi, and redshift, z, derived from recent ob- 



servations of the Lya forest ( Kim et al. 200l| ) and from 
Lyman Limit systems. The QSOs emissivity has been de- 
rived from the recent fits of Boyle et al. (2000| ), while 
we have used the stellar population synthesis model of 
Bruzual fc Chariot (2001 ) and a star- formation history 
from UV observations for the galaxy emissivity. Due to 
the present uncertainties in models and observations, we 
have used three values for the fraction of ionising photons 
that can escape a galaxy interstellar medium, /esc= 0.05, 
0.1 and 0.4, as suggested by local and high-z UV observa- 
tions of galaxies, respectively. 

The contribution of galaxies to the UV background 
is found to be comparable or larger than that of QSOs. 
This is consistent with other determinations of the galac- 



tic contribution to the background (Giallongo et al. 1997; 
Devriendt et al. 1998|; |Shull et al. 1999|). Taking into ac- 



count (in a rough fo rm) the contribution of rcemission 



from the IGM clouds ([Haardt fc Madau 1996| ; [Fardal et al 
199S), we find that all models with < /esc 0.10 can 



provide an ionising background within the limits measured 
from observations of the proximity effect at z 3 and 



softening of the UV background at z > 3. 

To summarise our result, a galaxy-dominated back- 
ground with /osc ^0.1 is consistent with the estimates of 
J(912A). In the hypothesis that the formation of struc- 
tures plays a neghgible effect, /esc ^ 0.05 — 0.1 is needed 
in order to explain the observed dN/dz of Lya absorbers. 
The rapidly improving knowledge derivable from numeri- 
cal simulations, the determination of the cosmological evo- 
lution of the Lyman forest, the proximity effect estimates 
of the ionising background, and the evolution of the in- 
tensity ratios of metal absorption lines will put soon con- 
straints on the relative galaxy/QSO contribution to the 
UV background and, together with direct measurements 
of the Lyman continuum emission from galaxies, will make 
it possible to address the issue of the evolution of the es- 
caping fraction of photons from galaxies as a function of 
z. 
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